Scaling of the Renyi entropies in gapped quantum spin systems: 
Entanglement-driven order beyond symmetry breaking 
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We investigate the scaling ol the Renyi a-entropies in one-dimensional gapped quantum spin 
models. We show that the block entropies with a > 2 violate the area law monotonicity and exhibit 
damped oscillations. Depending on the existence of a factorized ground state, the oscillatory behav- 
ior occurs either below factorization or it extends indefinitely. The anomalous scaling corresponds 
to an entanglement-driven order that is independent of ground-state degeneracy and is revealed by 
a nonlocal order parameter defined as the sum of the single-copy entanglement over all blocks. 
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During the past decade the use of entanglement in the 
study of complex quantum systems has developed at an 
increasingly fast pace [J-Q- Extended analyses have es- 
tablished the monotonic scaling of the von Neumann en- 
tropy in the ground state of one-dimensional spin mod- 
els 0, Q , and its profound relations with conformal field 
theory (CFT) in the gapless cases Much attention 

is currently being devoted to the Renyi entropies (RE), 
whose infinite hierarchy provides the most complete in- 
formation on the spectrum of the reduced density matrix. 
The general characterization of quantum states provided 
by the RE plays a relevant role in determining the scal- 
ing properties of numerical algorithms based on matrix 
product states UH)]. The RE are also a useful tool 



to determine the continuous or discontinuous nature of 
a phase transition [l4j and to estimate quasi-long-range 
order in low-dimensional systems 15[. Furthermore, the 
concept of topological entanglement entropy 3, 17 1 can 
be extended to the RE for which it has been shown to 
coincide with the total quantum dimension [l8| . 

Above all, there is a growing awareness that the entan- 
glement properties of quantum ground states provide the 
most fundamental characterization of quantum phases of 
matter beyond the traditional approach, borrowed from 
classical statistical mechanics, based on symmetry break- 
ing and local order parameters 0. In this perspec- 
tive, the RE and their topological components are be- 
ing actively investigated in various problems at the cut- 
ting edge of condensed matter physics, including Bosc- 
Hubbard spin liqui ds I20| , frustrated models on nontrivial 
lattice geometries (21| . non-Abelian fractional Hall sys- 



tems |22j , and low-dimensional gapless models [23J, \2A 



For a bipartite quantum system in a global pure state 
Pab =| ^ab){^ab I, the RE of subsystem A (and anal- 
ogously for subsystem B) arc defined as: 



S a (A) 



1 



ln[Tr(^)] , 



(1) 



where pa = T^bPab is the reduced density matrix of 
subsystem A, and the parameter a takes nonnegative real 
values. The von Neumann entropy Si is recovered in the 
limit a — > 1 + . In the opposite limit, the entanglement 
entropy = -ln(A mQX ), where X max is the largest 
eigenvalue of pa, and it coincides with the single-copy 
entanglement [25[. When the RE S a (A) are computed 
for blocks Aof£ spins in the ground state of gapless mod- 
els, it is observed that, for a > 2, they are characterized 
by large sub-leading corrections that violate monotonic- 
ity in £ [26rl28j. Investigations on the nature of these 
corrections have unveiled the role played by the parity 
of the number of spins forming the block [29143 1| . The 
significance of these findings stimulates the investigation 
of the scaling behavior of RE in general gapped systems. 

In this work we present a comprehensive study of the 
ground-state RE for one-dimensional gapped quantum 
spin models in external field h. We show that for models 
admitting ground-state factorization there exist two dif- 
ferent scaling regimes, separated by the factorizing field 
hf . In the region h < hf all the RE of order a > 2 exhibit 
a non monotonic scaling with damped oscillations in I. 
These oscillations are due to a series of crossovers be- 
tween the eigenvectors of the reduced block density ma- 
trix associated to its largest eigenvalue, and thus to the 
bipartite ground-state geometric block entanglement. In 
all models admitting factorization, the area law behavior 
of all RE is restored for h > hf. For models that do not 
admit a factorization point, the oscillatory behavior per- 
sists indefinitely, even beyond the critical point. These 
phenomena are associated to an entanglement-induced 
order revealed by a global order parameter defined as 
the sum of the single-copy entanglement over all block 
sizes. This entanglement-driven ordered pattern is in- 
dependent of the Ginzburg-Landau magnetic order and 
superimposed on it in a symmetry broken phase. 

To set the stage, let us first consider the class of trans- 



FIG. 1: RE S a as functions of the transverse field h for dif- 
ferent block sizes I in a bipartite, infinite XY chain, with 
anisotropy 7 = 0.4. Dotted black line: £ = 1; dot-dashed red 
line: £ = 2; solid green line: £ = 3; dashed blue line: i = 4. 
From the top left, clockwise, graphs of Si, S2, Ss, and Sao. 



lationally invariant, one-dimensional XY spin models 
ff,y = ^E(l+7)^^ + i+(l-7>W +1 -/iE< ' ( 2 ) 

i i 

where af (a — x,y,z) stands for the spin- 1/2 Pauli op- 
erator on site i, h is the external transverse field, and 
7 is the anisotropy, taking values in the interval [0,1], 
whose extremes correspond, respectively, to the gaplcss 
isotropic XX model and to the Ising case. In the thermo- 
dynamic limit, at h c = 1 the system undergoes a quan- 
tum phase transition, developing a nonvanishing spon- 
taneous magnetization M x . In the ordered phase the 
ground state | G) is twofold degenerate and the energy 
spectrum is gapped. For the XX model the spectrum 
is gapless for all < h < 1 and | G) is always non- 
degenerate. For every value of 7 the ground state | G) 
becomes f ully fac torized, i.e. product of single-site states, 
at hf = yl — 7 2 32H34| . The factorizing field hf is thus 
always below the critical field h c = 1 and coincides with 



it only in the isotropic and gapless case of the XX model. 

In Fig. Q]we report the RE S a (£, h) of parity preserv- 
ing ground states as functions of h for different values of 
£. The von Neumann entropy S\(£,h) is monotonically 
decreasing both in h and in £. For 1 < a < 2 and h < hf 
the RE violate monotonicity in h but remain monotoni- 
cally non-decreasing in I: S a (£+l,h) > S a (£,h). Finally, 
for a > 2 and h < hf all RE oscillate strongly in h, so 
that the monotonic scaling S a (£+ 1, h) > S a (£, h) is vio- 
lated. The number and amplitude of the oscillations be- 
come more pronounced, respectively, as £ and a increase. 
The oscillations are due to a series of crossovers between 
the eigenvectors corresponding to the largest eigenvalue 
Xmax of the reduced block density matrix. The number 
of crossovers is an increasing function of £, while the am- 
plitude is an increasing function of a. This phenomenon 
is due to the increasing weight of X ma x for increasing a. 



FIG. 2: Scaling of the RE in £ for different a and h, with 
7 = 0.4. Squares: a = 1; circles: a = 2; up-triangles: a = 8; 
down-triangles: a = 00. The first three panels (clockwise 
from top left) show the non monotonic behavior of S a (£, h) 
for q > 2 and three different values of h below the factor- 
ization point hf ~ 0.916. The fourth panel illustrates the 
restored monotonic behavior at h = 0.95 > hf, i.e. still be- 
low criticality but above factorization. 



In particular, for 1 = 1 the maximum of Soo is a parity 
inversion point at which the eigenvectors corresponding 
to Xmax swap between two Bell states of opposite parity. 

In Fig.[5]we report the RE as functions of £ for different 
values of a and h. The violation of the area law behavior 
for h < hf and a > 2 is characterized by damped oscil- 
lations in £. The implications of the anomalous scalings 
reported in Fig. [2] can be understood by recalling that 

the bipartite single-copy entanglement ^oo is monotonic 

(2) 

in the bipartite geometric entanglement £q ' = 1 — \ max 
[35J . Given a system of total size N, is the minimum 
distance between the ground state | G) and the set of 
pure bi-separable states | ipt) ®| 0jv-«), where | ipi) is a 
state of a block £, and | 4>N-e) is a state of the remain- 
der of the chain. From the last panel of Fig. [5] one sees 

(2) 

that above factorization £ G increases monotonically in 



t. On the contrary, when h < hf, Fig. [2] shows that 

£ G and all the RE of order a > 2 oscillate as a func- 
tion of £, with the amplitude of the damped oscillations 
increasing as h decreases. In particular, at sufficiently 
small fields the true ground state | G) is closest to the 
bi-separable state | V2) ® | 4>n-2)- This fact, together 
with the singlet-triplet inversion of the eigenvectors as- 
sociated to X m ax, implies that below factorization | G) 
tends to order in two-spin domains. The existence of 
this quasi-dimerization can be explained by looking at 
the behavior of the RE in systems with fully dimerized 
ground states such as the Majumdar- Ghosh chain (3f3 | 
and models with long-distance entanglement 3^, 38 1. In 
these chains, if the block boundary cuts a dimer, all the 
resulting block RE are strongly enhanced; viceversa, if 
the boundary falls between two different dimers, they 
are strongly suppressed. By increasing the block size £ 



3 



the two situations alternate themselves, giving rise to a 
permanent, undamped, oscillatory behavior of the block 
RE (including the von Neumann entropy) as functions of 
I. Comparing with gapped spin models in zero field (see 
Fig. [5J first panel) we observe an important analogy, i.e. 
the oscillatory behavior and the ensuing tendency to an 
arrangement in two-spin domains. However, the oscilla- 
tory behavior in gapped models holds only for the RE of 
order a > 2 and therefore the tendency to dimerization 
is only partial, as it involves exclusively the eigenstates 
of the block reduced density matrix associated to X ma x ■ 
For finite h (sec Fig. [2j second and third panels), the 
frequency of the oscillations decreases and, crossing fac- 
torization, the onset of a unique domain comprising the 
entire system restores the area- law monotonic scaling (sec 
Fig- [H fourth panel). Exactly at factorization (h = hf) 
the system undergoes a further phase change beyond the 
quantum phase transition occurring at the critical point 
h c ; in the region h < hf < h c an entanglement-driven 
order is established that is superimposed on the global 
magnetic one. As the oscillatory behavior is maximized 
by the single-copy entanglement Soo (£, h) , this global or- 
der due solely to entanglement is naturally character- 
ized by the nonlocal order parameter defined by summing 
Soo(£, h) over all block sizes: 



r = — minfS'o 



(£,h)~S oo (£-l,h),0}. (3) 



In Fig. |3] we report the contour plot of T for the class 
of XY models. It vanishes identically for h > hf. Ap- 
proaching the Ising limit it is strongly reduced, corre- 
sponding to strongly damped oscillations. Approaching 
the XX limit T increases and diverges exactly at 7 = 0, 
consistently with the fact that in gapless models below 
the critical point the oscillations persist indefinitely at 
any I. Eq. ([3]) and Fig. [3] show that such nonlocal order 
below factorization is exclusively due to the entanglement 
properties of each individual ground state and therefore, 
unlike magnetic order, it is not a consequence of ground- 
state degeneracy and symmetry breaking. 

We will now show that the entanglement-induced order 
extends to the general class of gapped XYZ models: 



H x 



2 ^ 

i.i 



■JyOiOf- 



-J^crf-h^a*. (4) 



Here J M are the spin-spin couplings along the fj, = x,y, z 
directions and, without loss of generality, we set J x = 
1 > I J y I , I J z I . These systems undergo a quantum phase 
transition at h = h c . The relevant difference between the 
general XYZ models and the XY chain lies in the fact 
that while the ground state of the latter is always sub- 
ject to a factorization at h = hf, the models described by 
Eq. (0| admit a factorized ground state at a finite value 
h = hf = J (J x + J z )(J y + J z ) if and only if J z > —J y 




FIG. 3: Nonlocal order parameter F as a function of h and 
7. Left inset: one-dimensional projection as a function of h 
for different values of 7. Red dotted line: 7 = 0.2; green solid 
line: 7 = 0.4; black dashed line: 7 = 0.65. Right inset: one- 
dimensional projection ELS Si cl function of 7 for different values 
of h. Red dotted line: h — 0.2; green solid line: h = 0.65; 
black dashed line: h = 0.87. In both insets the vertical grid 
lines correspond to the factorization points. 
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FIG. 4: Single-copy entanglement Sao as a function of h for 
different values of I in a XY Z model with J x = 1, J y = 0.7, 
J z = 0.3, and hf ~ 1.14. Squares: 1 = 1; circles: £ = 2; 
up-triangles: t = 3; down-triangles: I = 4. Upper right inset: 
von Neumann entropy Si. Lower left inset: Soo (circles) and 
Si (squares) as functions of £ at h = 0. 



H, 33|. The models of the class Eq. (0| are not exactly 
solvable. In order to determine the block reduced den- 
sity matrices we diagonalizcd the system by means of the 
Density Matrix Rcnormalization Group (DMRG) [^lifjj 
applied to open chains of up to 128 spins. For mod- 
els admitting a factorization point, the results of the 
XY case carry over essentially unmodified to the gen- 
eral XYZ instance. In Fig. [4] we report the scaling of 
the RE in XYZ models that possess a factorizing field 
hf (J z > ~Jy)- The qualitative behavior is analogous 
to that of the XY models and confirms the existence of 
an entanglement-induced order of two-spin domains for 
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FIG. 5: Single-copy entanglement Soo as a function of h for 
different values of t in a XYZ model with J x = 1, J y = —0.7, 
J, = —0.3, and no factorization point. Squares: I = 1; circles: 
1 = 2; up-triangles: £ = 3; down-triangles: t = 4. Upper 
right inset: von Neumann entropy Si. Upper left inset: S^o 
(circles) and Si (squares) as functions of £ at h = 0. 



fields h < hf also in the general case. The investigation 
of models that do not admit a factorized ground state 
at any finite value of the magnetic field is summarized 
in Fig. [5] At variance with the case admitting ground 
state factorization, the area-law order between blocks of 
different size S a (£ + l,h) > S a (£, h) appears to be always 
violated for all values of h, as one can see comparing 
the main plots of Figs. 0] and [5] On the other hand, 
a rigorous confirmation of this behavior for all values of 
h is beyond the current numerical possibilities, since the 
differences between RE of the same order a and differ- 
ent block size £ fall off extremely rapidly with h. In the 
absence of factorization all the RE behave smoothly and 
do not acquire local maxima. Moreover, considered as 
functions of the block size £ (at fixed external field h) 
they exhibit the same oscillating behavior as in models 
with factorized ground states (See the left inset of Fig. [5]). 
These features can be intuitively understood by consid- 
ering that in these models factorization, so to speak, is 
moved towards infinitely large values of h, corresponding 
to a classical saturation. This behavior confirms that the 
entanglement-induced order can exist also in a symmetry 
unbroken phase with vanishing local order parameters. 

In conclusion, we have investigated the Renyi entangle- 
ment entropies in one-dimensional gapped quantum spin 
models. We have showed that a violation of the area law 
scaling behavior, analogous to that occurring in the gap- 
less cases, holds for all Renyi entropies of order a > 2 
and for external fields h < hf < h c for models admitting 
ground-state factorization. For models that do not admit 
factorization, it appears to extend to all values of h. This 
anomalous scaling characterizes the existence of a novel 
type of entanglement-driven order. Unlike the magnetic 
order, this ordering does not rely on the onset of ground- 
state degeneracy, but is rather due to the entanglement 
properties of each parity preserving ground state. The 
factorizing field plays a key role in the understanding of 



the anomalous scaling of the RE, by clarifying the origin 
of this phenomenon even in the gaplcss case in which the 
factorization and critical points coincide. It also explains 
why the anomalous scaling is not observed in systems, 
like the Ising model, for which hf = 0. The maximal 
violation of the area law occurs in the case of the RE of 
order a = oo, the single-copy entanglement. This fact 
allows to introduce a nonlocal order parameter, defined 
as the sum of the single-copy entanglement over all block 
sizes, that is nonvanishing for all fields h < hf and van- 
ishes identically for h > hf. On the other hand, the 
single-copy entanglement is monotonic in the bipartite 

(2) 

geometric entanglement £ G , which has been shown to 
be a universal lower bound to ground-state frustration 
41| . This correspondence suggests the existence of an 
intimate relation between entanglement scaling and frus- 
tration of purely quantum origin that will be the subject 
of further analysis. In such a perspective, the present in- 
vestigation might be fruitfully extended to geometrically 
frustrated models and systems with topological order and 
fractional excitations. 
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